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Abstract 

Measurements of CP-violating observables in neutrino oscillation experi- 
ments have been studied in the literature as a way to determine the CP- 
violating phase in the mixing matrix for leptons. Here we show that such 
observables also probe new neutrino interactions in the production or detec- 
tion processes. Genuine CP violation and fake CP violation due to matter 
effects are sensitive to the imaginary and real parts of new couplings. The 
dependence of the CP asymmetry on source-detector distance is different 
from the standard one and, in particular, enhanced at short distances. We 
estimate that future neutrino factories will be able to probe in this way new 
interactions that are up to four orders of magnitude weaker than the weak 
interactions. We discuss the possible implications for models of new physics. 
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I. NEW CP VIOLATION IN NEUTRINO INTERACTIONS 



In the future, neutrino oscillation experiments will search for CP-violating effects PH2"2"fl . 
The Standard Model, extended to include masses for the light, active neutrinos, predicts 
that CP is violated in neutrino oscillations through a single phase in the mixing matrix for 
leptons. This effect is suppressed by small mixing angles and small mass differences. 

It is not unlikely, however, that the high-energy physics that is responsible for neutrino 
masses and mixing involves also new neutrino interactions. Such interactions provide new 
sources of CP violation. In this work we study CP-violating effects due to contributions 
from new neutrino interactions to the production and/or detection processes in neutrino 
oscillation experiments. We investigate the following questions: 

(i) How would new, CP-violating neutrino interactions manifest themselves in neutrino 
oscillations? 

(ii) Are the effects qualitatively different from the Standard Models ones? In particular, 
can we use the time (or, equivalently, distance) dependence of the transition probability to 
distinguish between Standard Model and new CP violation? 

(hi) How large can the effects be? In particular, do the new interactions suffer from 
suppression factors related to mixing angles and mass differences? 

(iv) Can the new CP violation be observed in proposed experiments? What would be 
the optimal setting for these observations? 

(v) Which models of New Physics can be probed in this way? 

The plan of this paper goes as follows. In section II we present a parameterization of 
the New Physics effects that are of interest to us and explain the counting of independent 
CP-violating phases in our framework. In section III we evaluate the New Physics effects on 
the transition probability in neutrino vacuum oscillation experiments. (A full expression for 
the transition probability, without any approximations concerning mixing angles and mass 
differences, is given in Appendix A.) In section IV we investigate the resulting CP asymmetry 
and compare the New Physics contribution to the standard one (that is, the contribution 
to the asymmetry from lepton mixing). In sections V and VI we evaluate the New Physics 
effects on, respectively, the transition probability and CP asymmetry, in neutrino matter 
oscillations. In section VII we study how these effects can be observed in future neutrino 
factory experiments. In particular, we estimate a lower bound on the strength of the new 
interactions that can be observed in these experiments. This lower bound is compared to 
existing model-independent upper bounds in section VIII. We summarize our results and 
discuss some of the implications that would arise if a signal is experimentally observed in 
section IX. 

II. NOTATIONS AND FORMALISM 

In this section we give a model-independent parameterization of New Physics effects 
on production and detection processes in neutrino oscillation experiments. We put special 
emphasis on CP-violating phases. 

We denote by i = 1,2,3, the three neutrino mass eigenstates. We denote by \u a ) 
the weak interaction partners of the charged lepton mass eigenstates a~ (a = e,/i, r): 
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\v a ) = ^2U ai \ui). 



(2.1) 



Whenever we use an explicit parameterization of the lepton mixing matrix p3|,p4|, we 



will use the most conventional one: 



u = u 23 u 13 u 12 
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sin 9ij and Cjj = cos 6^-. Alternatively, a convention-independent definition of the 



phase 5 that we will use in our calculations is given by 
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We consider new, possibly CP-violating, physics in the production and/or detection pro- 
cess. Such effects were previously studied in Ref . |^5| , and we follow closely the formalism of 
that paper. Most of the analysis in f2~5 |, however, was carried out assuming CP conservation. 
We parameterize the New Physics interaction in the source and in the detector by two sets of 
effective four-fermion couplings, (G^p)^ and (G^ P ) al 3, where a,j3 = e,fi,r. Here (Gfjp)^ 
refers to processes in the source where a up is produced in conjunction with an incoming a~ 
or an outgoing a + charged lepton, while (G^p)a/3 refers to processes in the detector where 
an incoming vp produces an a~ charged lepton. While the SU{2)l gauge symmetry requires 
that the four-fermion couplings of the charged current weak interactions be proportional to 
GpSa/3, new interactions allow couplings with a 7^ (3. Phenomenological constraints imply 
that the new interaction is suppressed with respect to the weak interaction, 



NpJa/3 



NPJ a/3 



(2.4) 



For the sake of concreteness, we consider the production and detection processes that are 
relevant to neutrino factories. We therefore study an appearance experiment where neutrinos 
are produced in the process /i + — > e + iy a u fl and detected by the process v$d — > fi~u, and 
antineutrinos are produced and detected by the corresponding charge-conjugate processes. 
Our results can be modified to any other neutrino oscillation experiment in a straightforward 
way. The relevant couplings are then (G^ P ) e ^ and (G^p)^- It is convenient to define small 
dimensionless quantities e^j in the following way: 
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NPje/3 



z ef3 
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Gf + (G^p) ee \ 2 + |((jNp)e^| 2 + \ {G : 



NP)er\ 



(G d 



Gf + (G^p) MAt | 2 + K^np)^! 2 + KGotOz-jtI 2 



(2.5) 



Since we assume that |e^'J < 1, we will only evaluate their effects to leading (linear) 
order. New flavor-conserving interactions affect neutrino oscillations only at (9(|e| 2 ) and 
will be neglected from here on. (More precisely, the leading effects from flavor-diagonal 
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couplings are proportional to e (flavor-diagonal) xe (flavor-changing) and can therefore be 
safely neglected.) 

We use an explicit parameterization for only two of the e's, with the following convention: 

<, = K\e iS % e% = \e%\e<. (2.6) 
Alternatively, we can define the phases S e and S' e in a convention-independent way: 

We would like to conclude this section with a comment on the number of independent 
CP-violating phases in our framework. It is well known that the three-generation mixing 
matrix for leptons depends, in the case of Majorana neutrinos, on three phases. Two of 
these, related to the fact that there is no freedom in redefining the phases of neutrino fields, 
do not affect neutrino oscillations and are therefore irrelevant to our discussion. The other 
one is analogous to the Kobayashi-Maskawa phase of the mixing matrix for quarks. The 
freedom of redefining the phases of charged lepton fields is fully used to reduce the number 
of relevant phases to one. Consequently, it is impossible to remove any phases from the e*'^ 
parameters. Each of these parameters introduces a new, independent CP-violating phase. 

For example, when we discuss v e — > oscillations, our results will depend on e* M , e^ e 
and the UeiU^ (i = 1, 2, 3) mixing parameters. This set of parameters depends on three 
independent phases, one of which is the 5 of Eq. ( |2.3| ), while the other two can be chosen 
to be 5 e and 5' e of Eq. fl2.7f) . This situation is illustrated in Fig. [I], where we show in 
the complex plane the unitarity triangle and the e s,d parameters that are most relevant to 
v e ~ > Vp oscillations. 



III. THE TRANSITION PROBABILITY IN VACUUM 

In this section we derive the expression for the transition probability in neutrino os- 
cillation experiments as a function of the mixing matrix parameters and the New Physics 
parameters. 

We denote by v s e the neutrino state that is produced in the source in conjunction with 
an e + , and by v% the neutrino state that is signalled by pT production in the detector: 



1) = E [ u * + + e s eT U Ti ] \u t ) 
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(3.1) 



(Note that the norm of the states so defined is one up to effects of (9(|e| 2 ), which we 
consistently neglect.) We obtain the following expression for the transition probability P efl = 
\{Vp\v%{t))\ 2 , where vl(t) is the time-evolved state that was purely v s e at time t = : 



e/j, 
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Our results will be given in terms of Am|, Ay and £y, which are defined as follows: 



m- 



m), 



Ay = Am 2 /(2£) 



AyL/2, 



(3.3) 



where -E 1 is the neutrino energy and L is the distance between the source and the detector. 

Equation (|3.2|) will be the starting point of our calculations. The full expression for 
P efl in vacuum is given in Appendix A and has been used for our numerical calculations 
described below. To understand the essential features of our analysis it is, however, more 
useful to do the following. First, we separate P etli into a Standard Model piece, P™, and a 

sAI ; ' p — 0). This is the contribution to 



New Physics piece, P^J- What we mean by P^ is P e ^ a p 



P efl from the Standard Model extended to include neutrino masses but no new interactions. 
In contrast, PfJ contains all the e ^-dependent terms. Second, since the atmospheric and 
reactor neutrino data imply that \U &3 \ is small and the solar neutrino data imply that 
Am 2 2 / Am 2 3 is small, we expand P^ to second order and P^jf to first order in \U e3 \ and 
Am\ 2 . 

For P^ we obtain: 



efi 



±x z 21 \U e2 \ l \U„ 2 



| 2 + 4sin 2 X3i|f/ e3 r|t/ M 3| 



+ Ax 2l sm2x 31 lle(U e2 U* 3 U; 2 U^) 
- 8x21 sin 2 x 31 lm(U e2 U* 3 U* 2 U fl3 ) 



(3.4) 



The first term is the well known transition probability in the two-generation case. The 
second term gives the well known transition probability in the approximation that Am\ 2 = 0. 
The last term is a manifestation of the Standard Model CP violation. 
For P e ^ p we obtain 
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(3.5) 



The last three terms in this expression are CP-violating and would be the basis for our 
results. 



IV. CP VIOLATION IN VACUUM OSCILLATIONS 

To measure CP violation, one will need to compare the transition probability P e[l eval- 
uated in the previous section to that of the CP-conjugate process, P^. The latter will 
be measured in oscillation experiments where antineutrinos are produced in the source in 
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conjunction with e~ and detected through /z + production. It is clear that a CP transforma- 
tion relates the production processes, //~ — > t~v a v a i and ji + — > e + z/ a z/ a /. As concerns the 
detection processes, Vpu — > yU + c/ and z/gd — > the situation is less straightforward. We 
have Gp^ oc (pfi~\p,~ u Vp d\vpn) and oc d Vp u\Vpp). The relation is through 

CP and crossing symmetry, but for a four-fermion interaction this is equivalent to a CP 
transformation. 

CP transformation of the Lagrangian takes the elements of the mixing matrix and the 
e-terms into their complex conjugates. It is then straightforward to obtain the transition 
probability for antineutrino oscillations. Our interest lies in the CP asymmetry, 

Acv = ^, (4.1) 

-m- 

where 



P± = P eix ±P m . (4.2) 

We quote below the leading contributions for 'short' distances, x%\ < 1. In some of the 
observables, we consider two limiting cases for \U e3 \: 

The "large" s 13 limit : 2:21/2*1 < \{U eZ U^) / {U e2 U^)\, 
The small s 13 limit : x 2 i/x 31 > ^U^U^)/ {U^U^- (4.3) 

The CP conserving rate P + is always dominated by the Standard Model. It is given by 

= {%xl 1 \U e zU^\ 2 large si 3 , . . 

+ {S^U^U^ 2 small s 13 . 1 ' ] 

The CP-violating difference between the transition probabilities within the Standard 
Model can be obtained from Eq. ( |3.4|) : 



P SM = -lex^lmiU^U^U^U^). (4.5) 

As is well known, CP violation within the Standard Model is suppressed by both the small 
mixing angle \U e3 \ and the small mass-squared difference Amf 2 . More generally, it is propor- 
tional to the Jarlskog measure of CP violation, J = Tm{U e 2U* 2 U^U^) . For short distances 
(x2i,xsx <C 1), the dependence of P^ M on the distance is L 3 . Since it is CP-violating, it 
should be odd in L. The absence of a term linear in L comes from the fact that the Stan- 
dard Model requires for CP to be violated, that all three mass-squared differences do not 
vanish, that is, P_ oc A2iA 31 A 3 2. In the limit £21/^31 <C {(U^U^) / (U&U^l, we obtain the 
following Standard Model asymmetry: 

^- 2 - im (S)- (4 - 6) 

In the small S13 limit, the standard CP violation is unobservably small. 

The CP-violating difference between the transition probabilities that arises from the 
New Physics interactions can be obtained from Eq. (|3.5|): 
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NP f -8x 3 iXm[[/*3^ 3 (e^ + e^)] large s 13 , , , 

- \ -8x 21 lm[U: 2 U^{e% + e^)] small s 13 . 1 J 

We learn that CP violation beyond the weak interactions requires only that either |Z7 e 3| or 
Amfj be different from zero, but not necessarily both. Also the dependence on the distance 
is different: for short distances, P NP oc L. From Eqs. Q4.4TJ and ( [4.7[ ) we obtain the following 
new physics contribution to the CP asymmetry: 

— X -Zm (f^) large s 13 , 
A^ P = { ^ ^ (4.8) 
-i-Jm if±j^ small s 13 . 



asai V U e2 U /l2 



The apparent divergence of A^p for small L is only due to the approximations that we 
used. Specifically, there is an (9(|e| 2 ) contribution to P + that is constant in L ||25|| , namely 
P + = 0(\e\ 2 ) for L — ► 0. In contrast, P_ = in the L — > limit to all orders in |e|. 
Equations (|4.7f) and (|4.8|) lead to several interesting conclusions: 

(i) It is possible that, in CP-violating observables, the New Physics contributions com- 
pete with or even dominate over the Standard Model ones in spite of the superweakness of 
the interactions (|e| <C 1). Given that for the proposed experiments x 3 i < 1, it is sufficient 



that 



max 



OUKJ) >min(|f/ e3 |,* 21 ), (4.9) 



for the new contribution to the CP-violating difference P_ to be larger than the standard 
one. 

(ii) The different distance dependence of P^ p and P_ will allow, in principle, an un- 
ambiguous distinction to be made between New Physics contributions of the type described 
here and the contribution from lepton mixing. 

(iii) The 1/L dependence of A^p suggests that the optimal baseline to observe CP vio- 
lation from New Physics is shorter than the one optimized for the Standard Model. 

We carried out a numerical calculation of the probabilities P± and asymmetry Aqp as a 
function of the distance between the source and the detector. We use E u = 20 GeV, which 
is the range of neutrino energy expected in neutrino factories. For the neutrino parameters, 
we take Amf^ = 3 x 10~ 3 eV 2 and tan 2 823 = 1, consistent with the atmospheric neutrino 
measurements [ 2lf| , and Am 2 x = 10~ 4 eV 2 and tan 2 9±2 = 1, consistent at present with the 



LMA solution of the solar neutrino problem [^,^]. As concerns the third mixing angle 
and CP-violating phase in the lepton mixing matrix, we consider two cases. First, we 
take S13 = 0.2, close to the upper bound from CHOOZ |28|j29|j26[1 , and 5 = ir/2. This set 



of parameters is the one that maximizes the standard CP asymmetry. Second, we take 
S13 = 0, in which case there is no standard CP violation in the lepton mixing. As concerns 
the effects of New Physics, we demonstrate them by taking only |e*| 7^ 0. With our first 
set of mixing parameters (maximal standard CP violation), we take \e s e \ = 10~ 3 and 5 e = 0. 
With our second set of mixing parameters (zero standard CP violation), we take \e s e \ = 10~ 4 
and 5 e = it/2. Our choice of CP-violating phases can be easily understood on the basis of 
Eq. ([4.8|) : in the large S13 limit, the CP asymmetry depends on a.ig[e s e ^/ (U e 3U* 3 )} = 5 e — 5, 
while in the small S13 limit it depends on a.ig[e s efi / (U e 2U* 2 )} = <5 e - We use the full expression 
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for the transition probabilities that is presented in Appendix A. Consequently, the only 
approximation that we make is that we omit effects of (9(|e| 2 ). 

The results of this calculation are presented in Fig. |2|. The left panels correspond to the 
first case (maximal standard CP violation) and the right ones to the second (zero standard 
CP violation). For each case we present, as a function of the distance between the source 
and the detector, P + (dotted line), P^ M and A^f> (dashed lines in, respectively, upper and 
lower panels), and P^ p and A^F P (solid lines in, respectively, upper and lower panels). 

We learn a few interesting facts: 

(i) The New Physics contribution to CP violation can dominate over even the maximal 
standard CP violation for values of |e| as small as 10~ 4 . This is particularly valid for distances 
shorter than 1000 km. 

(ii) The approximations that lead to Eqs. ( |3~? ) and ( |4.8| ) are good for L < 5000 km. 

(iii) As anticipated from our approximate expressions, for short enough distances, P^ p 
grows linearly with distances and A^, is strongly enhanced at short distances. 

(iv) In the large S13 case, the new CP violation is sensitive mainly to the phase difference 
5 — 5 € and is almost independent of the solar neutrino parameters. 

(v) In the very small S13 limit, the new CP violation is proportional to sin5 e . The rate 
P NP is suppressed by the solar neutrino mass difference and mixing angle. 



V. THE TRANSITION PROBABILITY IN MATTER 

Since long-baseline experiments involve the propagation of neutrinos in the matter of 
Earth, it is important to understand matter effects on our results. For our purposes, it is 
sufficient to study the case of constant matter density. Then the matter contribution to the 
effective v e mass, A = y/2GpN e , is constant. 

One obtains the transition probability in matter by replacing the mass-squared differences 
Aij and mixing angles U a i with their effective values in matter, A™ and U%£. The full 
expression for P efl in matter can then be written in terms of x™ and £7™ by a straightforward 
modification of the vacuum probability given in Appendix A. To understand the matter 
effects it is, however, more useful to take into account the smallness of \U e3 \ and x\2- We 
will expand the transition probability in these parameters to second order for Pjr 1 and to 
first order for P^ p . 

For the Standard Model case, we obtain: 




I cosx 31 TZe[U* 3 U^U e 2U* 2 } - smx 3 ilm[U* 3 U M3 U e2 U* 2 )\, 

where 

B = A 31 -A. (5.2) 

Again, the first term is the full result for two generations, and the second is the full result for 
the case of A 21 = 0. The last term violates CP. In the limit A = 0, Eq. (|3.4j ) is reproduced. 
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Note that our definition of B is such that B changes sign according to whether A31 is larger 
or smaller than A. This is different from the usual convention where B = |A 31 — A\. The 
Standard Model results are an even function of B and either definition can be used. But for 
the New Physics results given below, the choice of convention is important. 
For the New Physics contribution we find: 
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A31 
B 

A21 
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B 
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u: 3 u, 3 (eyi - \u, 



/u3 1 



v 2 y 

) Zm[f/; 2 ^ 2 (e^l^al 2 + ^^3^3)] }• 



+ 8sm I ~2~ J sm ( ) sm •''■•!] 
A31 

A 2 i 



(5.3) 



Unlike the case of vacuum oscillation, P_ will get contributions from both CP-violating 
terms (proportional to the imaginary parts of various combinations of parameters) as CP 
conserving terms (proportional to the real parts). 

Note that, in addition to the effects of new neutrino interactions in the source and in the 
detector, there could be other, independent effects due to new neutrino interactions with 
matter during their propagation PUH32]. Such effects have been studied in the context of 



solar and atmospheric neutrinos (see e.g. Refs. p3|-p5[]) but we neglect them here. 



VI. CP VIOLATION IN MATTER OSCILLATIONS 

Since matter in Earth is not CP symmetric, there will be contributions to Aqp even in 
the case when there is no CP violation. It is our purpose in this section to evaluate these 
contributions and, in particular, the fake asymmetry that is related to the real part of e. 
We denote the matter-related contribution to P by P m = P-(A) — P-(A = 0). Since the 
leading contributions to P + are the same as in the vacuum case [Eq. ( |4.4[ )], we can similarly 
define the matter-related contribution to Aq?\ Aq P = P™/P + . Note that in the evaluation 
of P S p, from the expressions that we found for P eM we need not only to replace U ai and e*'f 
with their complex conjugates, but also A with —A. 
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For the Standard Model, we obtain from Eq. ( |5.1|) , in the small X31 and large S13 limits, 

(P m ) SM = ™ x ii (£) \ u ^\ 2 - 

In the small s i3 limit (x 2 i/x 3 i ^> KU^U^) / (U^U^l) the Standard Model effect is unob- 
servably small, and we do not consider it here. Taking into account that [see Eq. ( [4.41) 1 

P+ « 8a4 1 E/esf/^a I 2 , we get 

(^cp) SM = ^31 ( • (6-2) 




3 31 VA 31 

For the New Physics contribution, we obtain from Eq. ( |5.3|) , in the small x 3 i limit, 

rp ^NP = (^i^MU^U^-elJ) larger, 

IS^iA^^^We-^)] small Sl3 , { - } 

and 

large s 13 , 

(6.4) 

small S13. 

We would like to make a few comments regarding our results here: 

(i) Each of the four contributions has a different dependence on the distance. In the 
short distance limit, we have 

(P m ) SM ocL 4 , P! M ocL 3 , (Pr) NP ocL 2 , P^ocL, (6.5) 

and, equivalently, 

(A™ P ) SM ocL 2 , A$g<xL, (A™ P ) NP oc L°, A^ocl/L. (6.6) 

One can then distinguish between the various contributions, at least in principle. 

(ii) If the phases of the e's are of order 1, then the genuine CP asymmetry will be larger 
(at short distances) than the fake one. 

(hi) It is interesting to note that the search for CP violation in neutrino oscillations will 
allow us to constrain both TZe(e) and Im(e). 

We carried out a numerical calculation of the probabilities P™ and asymmetry A™ P as a 
function of the distance between the source and the detector. We use again E u = 20 GeV, 
Am 2 ! = 3 x 10~ 3 eV 2 , tan 2 6 23 = 1, Am 2 21 = 10" 4 eV 2 , tan 2 12 = 1 and s 13 = 0.2 or 0. For 
the New Physics parameters, we take |e® | = 10 -3 . To isolate the matter effects we now, 
however, switch off all genuine CP violation, that is, we take 5 = 5 e = in both cases. 

The results of this calculation are presented in Fig. |3[ The left panels correspond to 
the first case (large Si 3 ) and the right ones to the second (vanishing si 3 ). For each case we 
present, as a function of the distance between the source and the detector, P + (dotted line), 
(P!") SM and (Aq P ) sm (dashed lines in, respectively, upper and lower panels), and (p^™) NP 
and (Aq P ) np (solid lines in, respectively, upper and lower panels). 
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We learn a few interesting facts: 

(i) The New Physics contribution to the fake CP violation can dominate over the standard 
contribution for values of |e| as small as 10 -4 . This is particularly valid for distances shorter 
than 500 km. 

(ii) As anticipated from our approximate expressions, for short enough distances (p™) NP 
grows quadratically with distances and (t4™ p ) np is independent of the distance. 

(hi) Both the standard and the new contribution to P™ are suppressed by a small S13. 
The S13 suppression is however stronger for P + than it is for (p™) NP . Consequently, the 
New Physics contribution to (t4™ p ) np becomes very large for vanishing s 13 . 

In reality, the measured P_ and Acp will be affected by both genuine CP-violating 
contributions and matter-induced contributions. This situation is illustrated in Fig. ||. We 
present P + (dotted curve), P^ M and (dashed curves in, respectively, upper and lower 
panels), and P^ p and A^p (solid curves in, respectively, upper and lower panels), as a 
function of the distance. For the neutrino parameters, we always take Am^ = 3 x 10~ 3 eV 2 
and tan ^23 — 1> consistent with the atmospheric neutrino data. For the other parameters, 
we take three cases: (a) Left panel: we take the LMA parameters (Ar^j = 10~ 4 eV 2 and 
tan #12 = 1), 'large' S13 = 0.2 and maximal phase 5 = tt/2. This choice of parameters gives 
maximal standard CP violation. For the New Physics parameters we take le^J = 10~ 3 and 
5 e = 0. (The reason for the choice of phase is that the dominant contributions depend on 
5 - 5 e .) (b) Middle panel: we take the SMA parameters (Am 2 x = 5.2 x 10~ 6 eV 2 H||27|, 
tan 2 0i2 = 7.5 x 10" 4 ), s 13 = 0.2, 5 = tt/2, \e s e A = 10~ 3 and 5 e = 0. Here the standard CP 
violation is unobservably small, but the standard matter effects are still large, (c) Right 
panel: we take the LMA parameters and Si 3 = 0. With a vanishing si 3 , the total transition 
probability is highly suppressed as is the standard matter effect, and standard CP violation 
vanishes. For the New Physics parameters we take je^J = 10 -4 and 5 t = tt/2. We take a 
smaller le^J so that our approximation will not break down. 

We would like to emphasize the following points: 

(i) Similar three cases will be the basis, in the next section, for our analysis of the 
sensitivity of CP-violating observables measured in neutrino factories to New Physics effects 
(see Fig. ||). 

(ii) With large Si 3 , the dependence of the New Physics effects (and of the standard 
matter-induced effects) on the solar neutrino parameters is very weak. 

(hi) A small or even vanishing S13 will suppress all the rates and will introduce a strong 
dependence on the solar neutrino parameters. The New Physics contributions to Aqp will 
be, however, only little affected because both the standard CP conserving rate and the New 
Physics CP-violating rate are suppressed in the same way. 

(iv) With large S13, the New Physics CP-violating effects are dominated by the combi- 
nation 6 — 5 e . With small (but not vanishing) S13, the dependence is on both 5 — 6 6 and 

Se- 

(v) For distances shorter than 800 km, the effects of |e| > 10~ 3 are always dominant. For 
distances shorter than 300 km, the New Physics dominates even for |e| ~ 10~ 4 . 
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VII. LONG BASELINE EXPERIMENTS 



We would like to quantify the sensitivity of a neutrino factory to the CP-violating effects 
from new neutrino interactions. For this purpose, we consider the measurement of the 
following integrated asymmetry 0: 

N\ji-]/N [e-]\ + -N\p+]/N Q [e+]\_ 



AcP N[f,-)/N [e-)\ + + N[f,+)/N [e+)\J {7A) 

Here iV[/i~]/iVo[e~]| + refers to an oscillation experiment that has decay as its produc- 
tion process: N[fi~] is the measured number of wrong-sign muons while iVo[e _ ] is the 
expected number of v e CC interaction events (in the absence of oscillations). Similarly, 
iV[/i + ]/A r [e + ]|_ refers to an oscillation experiment that has pT decay as its production pro- 
cess: iV[/i + ] is the measured number of wrong-sign muons while iVo[e + ] is the expected 
number of FjT CC interaction events (again, in the absence of oscillations). The measured 
number of wrong-sign muon events can be expressed as follows: 

N\pT]\ + = ^^ff / dEJ u (E u )a cc (E u )P e ,(E u ), (7.2) 

where Nt is the number of protons in the target detector, is the number of useful muon 
decays, is the muon energy and m M is the muon mass. The function f v {E u ) is the energy 
distribution of the produced neutrinos. We assume that the muons are not polarized, in 
which case f v {E v ) = 12a; 2 (l — x) with x = E u /E fl . Finally, acc(E u ) is the neutrino-nucleon 
interaction cross section which, in the interesting range of energies, can be taken to be 
proportional to the neutrino energy: o C c — °~oE u with cr = 0.67 x 10~ 38 cm 2 /GeV for 
neutrinos and <To = 0.34 x 10~ 38 cm 2 /GeV for antineutrinos. The expression for iV[e _ ]| + is 
obtained by an integral similar to Eq. (]7.2|) , except that P eA1 is replaced by 1. 

We define A^p as the contribution from new physics (that is, e-dependent) to the inte- 
grated CP asymmetry. We take into account both genuine CP-violating and matter-induced 
contributions. (In the limit of a real lepton mixing matrix, that is, no standard CP violation, 
the first contributions are proportional to Im(e) and the latter to lZe(e).) We define A A 
to be the statistical error on Aqp- In order to quantify the significance of the signal due to 
New Physics, we compute the ratio A^p/AA. 

The statistical error, AA, scales with distance and energy as follows: 

AA ~ 1 oc . 1 oc -=. (7.3) 

/iV[/i+]|_ + iV[ / i-]| + JPl M N cc ^ 



To find this scaling, we took into account that the number of CC interactions scales as 
Ncc oc El/L 2 while, for L < 3000 km, Pf M oc L 2 /E^. Consequently, the dependence of 
A A on the distance is very weak. Given our results for Aqp, we obtain the following scaling 
with distance of the signal-to-noise ratio: 

A^/AA oc I lj/L genuine CP-violating effects, ^ ^ 

c \ const (L) matter induced effects. 
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This behavior is illustrated in Fig. |5| where we display the signal-to- noise ratio, A^lp / A A, 
as a function of the distance. For simplicity, we consider only the effect of e s . The standard 
CP violation is presented only in the upper panel, where it corresponds to maximal A^f> 
(LMA parameters: Am^ = 10~ 4 eV 2 and tan #12 = 1, large S13 and S = vr/2), while the 
middle panel has unobservably small A^ (SMA parameters: Am^ = 5.2 x 1CT 6 eV 2 and 
tan 2 ^12 = 7.5 x 10~ 4 ), and the lower panel has zero A^p 1 (sx3 — 0). As concerns the new 
CP violation, the dashed line corresponds to the case with maximal CP-violating phase 
(5 e = |) and the solid line corresponds to purely matter-induced asymmetry (5 e = 0). In 
our calculations we have assumed a total of 10 21 useful \i~ decays with energy = 50 GeV 
and a 40 kt detector. 

It is clear from the figure that the maximal sensitivity to new, CP-violating contribu- 
tions to the production or detection processes will be achieved with shorter distances, while 
the sensitivity to CP conserving contributions through matter induced effects is almost 
independent of distance. 

A truly short baseline experiment can potentially probe the (9(|e| 2 ) CP conserving effects. 
But in this case, due to the small signal, systematic errors will dominate over the statistical 
ones discussed above. It is unlikely that |e| smaller than O(10~ 3 ) can be signalled in such a 
measurement. 

We next investigate the sensitivity to the size of the New Physics interaction that can be 
achieved by the measurement of the integrated CP asymmetry. In Fig. |6|, we show the regions 
in the \TZe{e s e ),Im(el )] plane that will lead to A^p/AA = 3 ( darker-shadow region) and 

App/Ai = 1 (lighter-shadow regions) at L = 732 km, the shorter baseline discussed for 
an oscillation experiment at a neutrino factory. We have assumed a total of 10 21 useful /i~ 
decays with energy = 50 GeV and a 40 kt detector. In all panels we have 5 = (no 
standard CP violation), Am^j = 3 x 10~ 3 eV 2 and tan 823 = 1, and the LMA parameters, 
Am^ = 10~ 4 eV 2 and tan #12 = 1. In the left panels we have S13 = 0.2 and in the right ones 
S13 = 0. In the upper panels Xm(e s efl ) > 0, which, for our choice of parameters, results in a 
constructive interference between the matter-induced and CP-violating effects, while in the 
lower panels Tm{e s e ^) < 0, which results in a destructive interference. 

In order to illustrate the expected improvement in sensitivity to the New Physics when 
the baseline is better optimized for this particular purpose, we plot in Fig. [7| the corre- 
sponding regions when the measurement of the integrated CP asymmetry is performed at a 
distance of L = 200 km. 

We would like to emphasize the following two points: 

(i) Fig. |] shows that |e| in the range 3 x 10~ 5 -10~ 4 would lead to a "3cr" effect. 

(ii) A shorter distance will improve the sensitivity to the new CP violation. Fig. [7] shows 
that, for 5 = 0, in which case CP-violating effects are proportional to Xm(e), an improvement 
by a factor of about 3 in the sensitivity to Im(e) is expected. In contrast, the sensitivity 
to lZe(e) is not affected by the choice of baseline since the new physics contribution to the 
matter-induced asymmetry is independent of L. 

(iii) A non-vanishing standard CP-violating phase, 5^0, together with a 'large' S13, will 
change the interference pattern between the matter-induced and CP-violating contributions 
from new physics. The reason is that now some of the contributions depend on S t — 5, so 
that TZe(e) and Im(e) do not correspond to matter- induced and CP-violating effects in any 
simple way. 
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VIII. PHENOMENOLOGICAL CONSTRAINTS 



The measurements of P e)1 and Pep are sensitive to the four effective couplings, e* , e s eTl 
e^ e and e^ T . These dimensionless couplings represent new flavor-changing (FC) neutrino 
interactions. They are subject to various phenomenological constraints. In this section, we 
present these bounds in order to compare them with the experimental sensitivity that we 
estimated in the previous section. 

The e s coupling gives the amplitude for the — ► e^u^u^ decay. For this process, 
there is no 5'[/(2) i -related tree-level decay that involves four charged leptons. Instead, by 
closing the neutrino lines into a loop, the four-Fermi coupling contributes to the fi — > ey 
and \i — > 3e decays. The question of how to extract reliable bounds from loop processes in 
an effective theory involves many subtleties. A calculation in the spirit of Ref. [Q yields 
very weak bounds. Instead, we quote here the bound in a specific full high energy model: if 
the effective iijJShV^^ coupling is induced by an intermediate scalar triplet, the constraint 
from the /i — > decay reads (see, for example, |3"7|] ) 

|ey < 5 x 10- 5 . (8.1) 

We emphasize again that the bound in (|8.1|) is model dependent and could be violated in 
models other than the one that we considered. 

The e s eT coupling gives the amplitude for the ^ — > e~h , T h , fl decay. The same coupling 
contributes also to the S77(2)i,-related process t~ — > jU + /i~e~. The experimental bound on 
the latter implies 

< 3.1 x 10~ 3 . (8.2) 

There could be SU(2) l breaking effects that would somewhat enhance the neutrino couplings 
with respect to the corresponding charged lepton couplings. These effects are discussed in 
detail in Refs. |38|j39|] where it is shown that they are constrained (by electroweak precision 
data) to be small. Since our purpose is only to get order-of-magnitude estimates of the 
bounds, we neglected the possible SU(2)l breaking effects in the derivation of (|8.2|) . 

The e^ e coupling gives the amplitude for u e d — > \i~u. It is constrained by muon conversion 

| e ; e |<2.1xl0- 6 . (8.3) 

The e^ T coupling gives the amplitude for v T d — > n~u. It is constrained by the r~ — > \T p 
decay [|9j: 

l<rl<10" 2 . (8.4) 

The bound on |e~ r | is the weakest that we obtain. Moreover, it is not unlikely that it is 
indeed the largest of the couplings since it is the only one not to involve a first-generation 
lepton. For precisely the same reason, however, its contribution to P eA1 is suppressed by 
an additional power of \U e ^\, which is the reason that it is omitted in our approximate 
expressions. 

Let us also mention that there is a generic bound of (9(0.1) on the purely leptonic 
couplings e s a a from universality in lepton decays and a somewhat weaker bound of (9(0.2) on 
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the semi-hadronic couplings e£g from universality in pion decays |39 |. While universality is 



experimentally confirmed to high accuracy, these bounds are rather weak because deviations 
from universality are 0(e 2 ). 

To summarize, we expect that all the e's that play a role in the transition probabilities 
of interest are of (9(1CT 3 ) or smaller. In the previous section, we learnt that proposed 
experiments might probe these couplings down to values as small as O(10~ 4 ). This means 
that the possibility to measure new neutrino interactions through CP violation in neutrino 
oscillation experiments is open. Conversely, such future experiments can improve the existing 
bounds on FC neutrino interactions which, at present, come from rare charged lepton decays. 



IX. CONCLUSIONS AND DISCUSSION 

We summarize the main points of our study: 

(i) CP-violating observables are particularly sensitive to new physics. The reason is 
that the standard CP violation that comes from the lepton mixing matrix gives effects that 
are particularly suppressed by small mass differences and mixing angles. Some of these 
suppression factors do not apply to new contributions. 

(ii) The fact that matter effects contribute to CP- violating observables means that these 
observables are sensitive to both the CP conserving and the CP- violating contributions from 
New Physics. 

(iii) The effects of New Physics in the production and detection processes depend on the 
source-detector distance in a way that is different from the standard one. One consequence 
of this situation is that, at least in principle, it is possible to disentangle standard and 
new effects. Another consequence is that in short distance experiments the new effects are 
enhanced. 

(iv) Our rough estimate is that future neutrino factories will be able to probe, through 
CP-violating observables, effects from new interactions that are up to about four orders of 
magnitude weaker than the weak interactions. 

(v) The sensitivity to New Physics effects is better than most of the existing model- 
independent bounds. 

We would like to mention that a similar (and, for specific models, even stronger) level 
of sensitivity may be achieved by other experiments that search for lepton flavor violation. 
Particularly promising are those involving muon decay and conversion (for a recenet review, 
see |J44| ): for example, a future experiment at PSI will be sensitive to B(fi — > ej) at the 
10~ 14 level ^5|, and the MECO collaboration has proposed an experiment to probe fi — e 



conversion down to 5 x 10 17 , four orders of magnitude beyond present sensitivities [HB 



If these experiments observe a signal, the search for related CP violation will become of 
particular importance. 

What type of new physics will be implied in case that a signal is observed? The e 
couplings represent effective four-fermion interactions coming from the exchange of heavy 
particles related to New Physics. If the New Physics takes place at some high scale Anp, 
then one can set an upper bound: 

(9-1) 

iv NP 
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The source of this bound is in the definition of e, which is the ratio of the four-fermion 
operator to Gp, and the fact that it is maximal when the New Physics contribution comes 
at tree level and the couplings are of order one. Since the expected experimental sensitivity 
is to |e| > O(10~ 4 ), we learn that we can probe models with 

A NP < 10 TeV. (9.2) 

If the New Physics contributes to the relevant processes only at the loop level, there is 
another suppression factor in |e| of order That would mean that such models can 

be probed only if Anp < 1 TeV. Finally, if the flavor changing nature of the interaction 
introduces a suppression factor, e.g. |e® J ~ m^/A^p, that by itself would be enough to make 
it unobservable in near future experiments. We thus learn that CP violation in neutrino 
oscillation experiments will explore models with a scale that is, at most, one to two orders of 
magnitude above the electroweak breaking scale, and where the flavor structure is different 
from the Standard Model. 

Another point concerns the Dirac structure of the four-Fermi interaction. We did not 
present it explicitly in our discussion of the G^p couplings. However, it is implicitly assumed 
in our discussion that the Dirac structure is the same as that of the weak interactions, i. e. 
a (V-A)(V-A) structure. The reason for that is that the effects that we discuss are a 
consequence of interference between weak and new interactions. A different Dirac structure 
would give strong suppression factors related to the the charged lepton masses. While 
our formalism would still apply, these suppression factors would make the related effects 
practically unobservable. 

We conclude that a signal is likely to imply new physics at a relatively low scale (up 
to 1-10 TeV) with new sources of flavor (and, perhaps, CP) violation. We know of sev- 
eral well motivated extensions of the Standard Model that can, in principle, induce large 
enough couplings. In particular, we have in mind loop contributions involving sleptons and 
gauginos in supersymmetric models, tree contributions involving charged singlet sleptons in 
supersymmetric models without R-parity, and tree contributions involving a triplet scalar 
in left-right symmetric models. In another class of relevant models, such as the model of 



Ref. ||40|| , active neutrinos mix with singlet neutrinos. (Here there can be Z-mediated contri- 



butions to the non-standard couplings, and the phenomenological constraints are different 



4l| , |42[| .) A detailed analysis of new neutrino interactions within relevant extensions of the 
Standard Model is beyond the scope of this paper, but preliminary results show that large 
enough couplings are allowed and in some cases even predicted C3| . 
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APPENDIX A: TRANSITION PROBABILITY IN VACUUM 



Neglecting terms of 0(e 2 ) and with no other approximations, we obtain the following 
expression for P eil : 

P e , = 4sin 2 :r 21 {|^ 2 | 2 |f/ e2 | 2 - 7fe[ el^U^U^ 2 + U^U^U^ 2 ) 

+ e%,{U; 2 U e2 \U el \ 2 + U^Ueil 2 ) + ^(U^UfaU^Un + U^U^U^) 

+ et T (u; 2 u e2 u* el u Tl + u^UaU^u*) - u: 2 u e3 u; 3 u, 2 



+ 2 sin 2x2i Itu 



+ 4e(u; 2 u e2 \u el \ 2 - ^mu^ 2 ) + et^u^u;^ - u^u^u;^) 
+ - u^u A uyj^) + u* e2 u e3 u; 3 u, 2 

+ Asin 2 x 31 {lM 2 |£/ e3 | 2 - Ke^U^U^ 2 + t/^al^il 2 ) 

+ e%,{U^U e3 \U el \ 2 + ^f/exl^l 2 ) + eUU^U^Urs + U^U^Un) 

+ e^u^Uesu:^ + u^u A uyj^) - u* e2 u e3 u; 3 u^ } 



+ 2 sin 2^31 Im 



+ <e(^3^3|f/el| 2 " U^U el \U e3 \ 2 ) + ^(U^U^U^ ~ U^U^U^) 

+ 4r( u U u e3 u :i u ri - UflU el uxJrt) - u* e2 u e3 u; 3 u, 2 



— 4 sin 2 x 32 Tie 



<,(u: 2 u,2\u, 3 \ 2 +u: 3 u, 3 \u, 2 \ 2 ) 



+ 4e(KMUe2\ 2 + u; 2 u e2 \u e3 \ 2 ) + e s eT (u* e2 u, 2 u; 3 u T3 + u* e3 u, 3 u; 2 u r2 ) 
+ ^Au; 3 u e3 u* e2 u T2 + u; 2 u e2 u: 3 u r3 ) + u* e2 u e3 u; 3 u, 2 



+ 2 sin 2x 32 Xm 



<M2U,2\u, 3 \ 2 - u* e3 u, 3 \u, 2 \ 2 ) 



+ 4e(u; 3 u e3 \u e2 \ 2 - u; 2 u e2 \u e3 \ 2 ) + e s eT (u* e2 u, 2 u; 3 u T3 - u: 3 u, 3 u; 2 u r2 ) 
+ el T {u; 3 u e3 u: 2 u T2 - u; 2 u e2 u: 3 u T3 ) + u* e2 u e3 u; 3 u, 2 . 



(Al) 
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FIG. 1. The neutrino parameters that dominate P e ^ in the complex plane. We 
show the relevant unitarity triangle, which is the geometrical presentation of the relation 
UelU*x + Ue2U* 2 + U^U*^ = 0, and the two parameters that describe the New Physics in the 
production, e® , and in the detector, e d * e . The three independent phases defined in the text, 5, 5 e 
and 5' e , are shown explicitly. The standard convention puts U e iU*i on the real axis. 
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L (Km) L (Km) 

FIG. 2. Transition probabilities and CP asymmetries in vacuum as a function of the distance. 
In the upper panels the curves correspond to P^M (dotted^ pSM (dashed) anc [ pNP ( so iid). In the 
lower panels the curves correspond to A^p (solid) and Afjp (dashed). In the left panels, S13 = 0.2, 
5 = it/2, \e s e J = 10~ 3 and 5 e = 0. In the right panels, S13 = 0, (e^J = 10~ 4 and 5 e = tt/2. In all 
curves E v = 20 GeV, Amj 3 = 3 x 10 -3 eV 2 , tan 2 9 23 = 1, Am^ = 10 -4 eV 2 and tan 2 6x2 = 1. 
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FIG. 3. Transition probabilities and fake CP asymmetries in matter as a function of the dis- 
tance. All CP-violating phases are set to zero. In the upper panels the curves correspond to 
pSM (dotted), (p™) SM (dashed) and (p™) NP (solid). In the lower panels the curves correspond 
to ( j 4q P ) np (solid) and (Aqp) sm (dashed). In the left panels S13 = 0.2, and in the right panels 
S13 = 0. In all curves E v = 20 GeV, Am| = 3 x 10~ 3 eV 2 , tan 2 23 = 1, Ara^ = 10" 4 eV 2 , 
tan 2 O12 = 1, 6 = 8 e = and \e° \ = 10" 3 . 
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FIG. 4. Transition probabilities and fake CP asymmetries in matter as a function of the 
distance. In the upper panels the curves correspond to P°M (dotted), P^ M (dashed) and 
(solid). In the lower panels the curves correspond to A^,p (solid) and A^p (dashed). In all curves 
E v = 20 GeV, Am| = 3 x 10" 3 eV 2 and tan 2 6> 23 = 1. In the left panels Ara^ = 10~ 4 eV 2 , 
tan 0i2 = 1, S13 = 0.2, S = ir/2, \e s efJ } = 10~ 3 and 5 t = 0. In the middle panels Am^ = 5.2 x 10" 6 
eV 2 , tan 2 6»i2 = 7.5 x 10 -4 , s 13 = 0.2, 5 = vr/2, |e^J = 10" 3 and 5 e = 0. In the right panels 
Am 2 .! = 10" 4 eV 2 , tan6>i 2 = 1, s 13 = 0, \e s \ = 10" 4 and 5 t = vr/2. 
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FIG. 5. The signal-to-noise ratio, ^4^p/A^4, as a function of the distance L. We considered the 
following parameters for the experiment: = 50 GeV, 10 21 decays and a 40 kt detector, and 
the following neutrino parameters: 5 = 0, Am| = 3 x 10 -3 eV 2 , tan #23 = 1. In the upper and 
lower (middle) panels we use the LMA (SMA) parameters. In the upper two (lower) panels we use 
S13 = 0.2(0). For the New Physics we take |e*| = 10~ 3 and 5 e = or tt/2. In the upper panel, the 
dotted curve gives the SM matter-subtracted asymmetry Aqp(5 = tt/2) — Aqp(S = 0). 
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FIG. 6. Regions in the plane of \TZe{e s e ^),Tm{e s e ^)] that give Aqp/AA = 3 (darker shadow) and 
1 (light shadow). For the experiment, we take L = 732 km, = 50 GeV, 10 21 fj,~ decays and a 
40 kt detector. For the neutrino parameters, we take 5 = 0, Am^ = 3 x 10 -3 eV 2 , tan 2 #23 = 1 ; 
Am|i = 10~ 4 eV 2 and tan 2 #12 = 1- In the left (right) panels we have S13 = 0.2(0). 
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FIG. 7. Regions in the plane of \TZe{e s eil ) ,Im(el )] that give A^p/AA = 3 (darker shadow) and 
1 (light shadow). For the experiment, we take L = 200 km, = 50 GeV, 10 21 fi~ decays and a 
40 kt detector. For the neutrino parameters, we take 5 = 0, Am^ = 3 x 10~ 3 eV 2 , tan 2 623 = 1, 
Am 2 ,! = 10~ 4 eV 2 and tan 2 Q\2 = 1- In the left (right) panels we have S13 = 0.2(0). Note that the 
scales in the right panels are different from the left panels and from Fig. |6|. 
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